Three Dimensional Relativistic Electromagnetic Sub-cycle Solitons 
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Three dimensional (3D) relativistic electromagnetic sub-cycle solitons were observed in 3D 
Particle-in-Cell simulations of an intense short laser pulse propagation in an underdense plasma. 
Their structure resembles that of an oscillating electric dipole with a poloidal electric field and a 
toroidal magnetic field that oscillate in-phase with the electron density with frequency below the 
Langmuir frequency. On the ion time scale the soliton undergoes a Coulomb explosion of its core, 
resulting in ion acceleration, and then evolves into a slowly expanding quasi-neutral cavity. 
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Time evolution of a three dimensional (3D) nonlinear 
wave, in general, differs drastically from that of ID or 
2D waves, as exemplified by the problems of the wave 
collapse jjj and of the transverse stability of solitons 
[pj. Relativistic electromagnetic solitons are now rou- 
tinely observed in 2D Particle-in-Cell (PIC) simulations 
H f|, ||, |(| in the wake of an intense short laser pulse 
propagating in an underdense plasma. Solitons attract a 
great attention because they are of fundamental impor- 
tance for nonlinear science Q| and are considered to be a 
basic component of turbulence in plasmas || . Thus the 
numerical identification of solitons, among the different 
kinds of coherent structures that are formed by an intense 
laser pulse in a plasma, stimulated a renewed interest in 
developing an analytical model |)| [l0| and in envisaging 
ways of detecting solitons experimentally ]ll[ . 

As was stressed in Ref. fL2| , a detailed description of 
the strong electromagnetic wave interaction with plasmas 
represents a formidable difficulty for analytical methods, 
due to the high dimensionality of the problem, the lack 
of symmetry and the importance of nonlinear and ki- 
netic effects. On the other hand, powerful methods for 
investigating the laser-plasma interaction have become 
available through the advent of modern supercomput- 
ers. In the case of an ultra-short relativistically strong 



plasm 3 



Soliton 



Laser pulse 



FIG. 1: Iso-surface of electromagnetic energy density corre- 
sponding to dimensionless value of (E 2 + B 2 )/Siv — 0.09/8-7T 
at t = 33.75 x 2-k/u. 



laser pulse, simulations with 3D PIC codes provide a 
unique opportunity for describing the nonlinear dynam- 
ics of laser plasmas adequately, including the generation 
of coherent nonlinear structures, such as the relativistic 
solitons. In this Letter we present numerical identifica- 
tion of a 3D subcycle relativistic soliton and complex 
spatial structure of its electromagnetic fields. 

Briefly summarizing the recent results in the analyti- 
cal and numerical investigation of relativistic solitons in 
plasmas we recall the development of the analytical the- 
ory of intense electromagnetic solitons [§, 10, 131, see also 
references in Ref. p3. The solitons found in ID and 2D 
simulations consist of slowly or non propagating electron 
density cavities inside which an electromagnetic field is 
trapped and oscillates coherently with a frequency be- 
low the unperturbed plasma frequency and with a spatial 
structure corresponding to half a cycle. One-dimensional 
sub-cycle relativistic electromagnetic solitons in an un- 
derdense plasma were observed for the first time in a 
PIC simulation in Ref. |^5[ , where the mechanism of 
soliton formation and the structure of the circularly po- 
larized soliton were investigated. The mechanism of soli- 
ton formation is related to the fact that the frequency 
in the rear part of an intense laser pulse propagating 
in an underdense plasma decreases down to the local 
Langmuir frequency because the pulse loses its energy 
while the number of photons is conserved. As a re- 
sult, the low-frequency part of the electromagnetic ra- 
diation of the pulse, propagating with very low velocity, 
is trapped inside the cavity in the electron density, and a 
sub-cycle soliton is formed. An exact analytical solution 
of the electron fluid - Maxwell equations representing ID 
circularly-polarized relativistic electromagnetic sub-cycle 
soliton was obtained in Ref. || in perfect agreement with 
ID PIC simulations. The 2D relativistic electromagnetic 
sub-cycle solitons discovered in the PIC simulations Q 
consist of two "pure" types of solitons: S'-solitons with 
transverse electric field and azimuthal magnetic field with 
respect to the symmetry axis, and P-solitons with the op- 
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FIG. 2: Soliton structure in the x-z and y-z planes at x = 14.5. In the upper row arrows represent electric field, E y and 
E z , where arrow length in the right of the figure corresponds to eE/m e ujc — 1, and background color indicates magnetic field 
eB x /m £ ujc. In the bottom row, electron density is shown with brightness scale, (a), (b) corresponds to t = 39.3 x 2n/u} and 
(c), (d) t = 40.2 x 2ty/lu. 







FIG. 3: Three dimensional structure of electric field (a),(b) 
and magnetic field (c),(d) in the soliton corresponding to Fig. 
2. In the electric field red and blue denote positive and neg- 
ative values of its divergence. In the magnetic field red and 
blue denote counterclockwise and clockwise rotation. Arrows' 
length corresponds to the magnitude of the fields. Frames 
(a) and (c) correspond to t = 39.3 x 2-k/ui; (b) and (d) to 
t = 40.2 x 2-k/u. 
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FIG. 4: Time dependence of the z-component of electric field, 
eE z /m e uc, at the center of the soliton, (a;, y, z) = (14.5, 0, 0). 




-2 

10 




b 


-3 

10 






4 

10 






( 


) 10 20 


30 








-2 

10 






-3 

10 






-4 







12 3 4 5 



FIG. 5: Ion acceleration during the soliton explosion: (a) the 
y-component of ion momentum, p y /m e c, and energy spec- 
trum of ions (b) in the domain [14.3:14.7,-2:2,-0.5:0.5] and 
(c) in the soliton core [14.3:14.7,-0.5:0.5,-0.5:0.5] at different 
normalized times uit. 
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FIG. 6: Three dimensional view of electron density at t — 
76x2tt/uj (a); and cross-sections of electron density in the y— z 
plane averaged over the space 35 < x < 37 at u)t/2n = 62.8 
(b), at 64.8 (c), and at 134 (d). 



posite structure: transverse magnetic field and azimuthal 
electric field. In contrast to electron vortices, which move 
across a density gradient, solitons move along the den- 
sity gradient towards the lower density. When a soliton 
reaches some critical density, it radiates its energy in the 
form of a low-frequency short electromagnetic burst [Q. 
The interaction of two 2D S-solitons leads to their merg- 
ing and the resulting soliton acquires the total energy of 
the two merged solitons Moreover, in an electron- 

ion plasma a 2D soliton evolves into a postsoliton j(| on 
an ion timescale due to the ion acceleration caused by the 
time-averaged electrostatic field inside the soliton. This 
effect leads to the formation of slowly expanding bubbles 
in the plasma density |(| . Note that here we use the term 
" soliton" for brevity, even if in principle the merging and 
the expansion on the ion time violate the strict definition 
of these structures as solitons. 

We present the results of a three dimensional simu- 
lation of laser induced sub-cycle relativistic electromag- 
netic soliton. We use REMP - Relativistic Electro- 
Magnetic Particle-mesh code based on the Particle-in- 
Cell method. This parallel and fully vectorized code 
exploits a new scheme of current assignment (lj that 
reduces unphysical numerical effects of the PIC method 
significantly. In the simulation the laser pulse propagates 
along the x-axis. The pulse is linearly polarized in the 



direction of the z-axis and its dimensionless amplitude 
is a = eE z / {m e u>c) — 1, corresponding to the peak in- 
tensity / = 1.38 • 10 18 W/cm 2 for the A = 1/mi laser. 
The laser pulse has a gaussian envelope with FWHM 
size 8A x 5A x 5A. Its focal plane is placed in front of 
the plasma slab at the distance of 3A. The length of the 
plasma slab is 13A. The plasma density is n e = 0.36n cr . 
Ions and electrons have the same absolute charge, and 
the mass ratio is mi/m e — 1836. The simulation box has 
660 x 400 x 400 grids with a mesh size of 0.05A. The 
total number of quasiparticles is 426 • 10 6 . The boundary 
conditions are periodic along the y- and z-axes and ab- 
sorbing along the x-axis for both the EM radiation and 
the quasiparticles. The simulations were performed on 
16 processors of the NEC SX-5 vector supercomputer in 
Cybermcdia Center, Osaka University. The simulations 
results are shown in figures where the space unit is 
the wavelength A of the incident radiation. 

In Fig. [j] we see one isolated soliton and a soliton 
train behind the laser pulse. A substantial part of the 
laser energy (up to 30%) is transformed into these co- 
herent entities. Figures. || and || show the structure of 
the isolated soliton with the electric and magnetic fields 
and the electron density at two different times with the 
interval approximately half of a soliton oscillation pe- 
riod. In the figures the space unit is the wavelength A 
of the incident laser pulse. The soliton consists of oscil- 
lating electrostatic and electromagnetic fields confined in 
a prolate cavity of the electron density. The cavity size 
is approximately 2A x 2A x 3A. The cavity is generated 
by the ponderomotive force and the resulting charge sep- 
aration induces a dipole electrostatic field. As seen in 
Fig. @, the charge density in the soliton oscillates up and 
down in the z-direction: at wt/2ir = 39.3 the electron 
hole is in the upper part of the figure and the elctron 
hump is in the lower, and vice versa - at totjlix = 40.2. 
The electric field at the soliton center is perpendicular 
to the direction of the laser propagation, so this mode 
differs from the laser-driven plasma wake. This field is 
so large that the quivering distance of electrons in the 
^-direction is of the order of the cavity size. This in 
turn results in continuous oscillations of the cavity. The 
soliton resembles an oscillating electric dipole. The oscil- 
lating toroidal magnetic field, shown in Fig. 0, indicates 
that besides the strong electrostatic field, the soliton also 
has the electromagnetic field. The electrostatic and elec- 
tromagnetic components in the soliton are of the same 
order of magnitude. 

Figure. || shows that the electric field in the soliton 
is poloidal, while the magnetic field is toroidal. We note 
that the magnetic field is mostly counterclockwise in the 
upper part of the soliton, and clockwise in the lower part. 
This structure of the electromagnetic field can be consid- 
ered as that of the lowest eigenmode of a cavity resonator 
with a deformable wall, and thus we call this structure a 
Transverse Magnetic (TM) soliton. Figure. ^ shows the 
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z-component of the normalized electric field at the center 
of the soliton. The electromagnetic field trapped in the 
oscillating cavity pulsates at the same frequency fig that 
is smaller than the surrounding unperturbed Langmuir 
frequency, w 0.87ui pe . The density of the cavity walls 
Therefore the electromagnetic energy can not 
be radiated away and, in addition, the soliton oscillation 
does not resonate with plasma waves. 

In the equatorial plane the structure of the three di- 
mensional soliton is similar to that of a two dimensional 
S-soliton, while that in the perpendicular planes is simi- 
lar to a two dimensional P-soliton. Considering the soli- 
ton as a wave packet, it has only half of one cycle in 
space, so we use the term "sub-cycle soliton". The dy- 
namic of the 3D soliton is clearly seen in the animations 
produced from the data (420 stills with period 0.05, see 
authors' website). Although as shown in Fig. 4 the field 
amplitude inside the soliton decreases because of its en- 
ergy losses due to ion acceleration and to the digging of a 
hole in the ion density as discussed below, the soliton life- 
time is sufficiently long to distinguish it from the other 
nonlinear modes generated by the laser pulse, such as the 
pulse wakefield or vortices. 

On the ion time-scale the soliton evolves into a postsoli- 
ton H . The ponderomotive force displaces the electrons 
outward and the Coulomb repulsion in the electrically 
non neutral ion core pushes the ions away with a pro- 
cess similar to the Coulomb explosion. The evolution 
of the wave-plasma interaction discussed can also be in- 
terpreted as a phenomenon similar to the wave collapse 
that is however saturated because the electrons are al- 
most completely evacuated by the strong wave. As the 
soliton amplitude decreases, the ions acquire a radial mo- 
mentum, as shown in Fig. |^. In contrast to the 2D S- 
polarized soliton discussed in Rcf. [||, the explosion of 
the 3D TM soliton is strongly anisotropic (the postsoli- 
ton cavity is elongated in the z-direction) and the energy 
spectrum of the accelerated ions has a minimum at zero. 
In Fig. | we see also the ion implosion near the center 
of the post-soliton. The ion implosion may appear in the 
post-soliton regime due to the electron heating by the 
trapped electromagnetic field at the density cavity walls, 
which causes the plasma ablation towards the cavity cen- 
ter similarly to the ion implosion and the dense plasma 
filament formation at the axis of the self-focusing channel 
discussed in Ref. |]l8) . 

At the last stage of the soliton evolution we see a slowly 
expanding postsoliton where the walls of the plasma cav- 
ity move with velocity v ~ 3T0~ 3 c. We notice that in the 
case of immobile ions we also see the soliton formation 
but the lifetime of the soliton is significantly longer and is 
determined by energy conversion into fast electrons sim- 
ilar to Landau damping. The almost isolated structures 
in the soliton train, which we consider as solitons even if 
they are not properly separated, tend to merge and form 
a foam of bubbles with relatively high-density (« 3n cr ) 



walls. 

A further proof of the electromagnetic nature of the 
solitary structure discussed above is provided by the 
3D PIC simulation in an inhomogeneous plasma with 
the density gradient in the w-direction from n e = to 
0.168n cr , as shown in Fig. We show that the soli- 
ton can propagate as a whole due to its electromagnetic 
nature, in contrast to the wake-field that remains at the 
same place due to its zero group velocity. In this case 
the dimensionless amplitude of the incident laser pulse 
is a = 3 and its FWHM size is 5A x 8A x 8A. The size 
of the plasma slab is 45A x 43A x 32A, and that of the 
simulation box is 90A x 48A x 32A. Ions are immobile. 
Initially, the laser pulse propagates along the x-axis, its 
symmetry axis intersects the plasma- vacuum interface at 
(x,y,z) = (5,26,0), where the local Lagmuir frequency 
is Upe — 0.3w. In Fig. [j] (a), we see the wakefield and 
a well-pronounced solitary structure. Frames (b) and (c) 
show half a period of the soliton evolution, similarly to 
Fig. H Frame (d) shows that the solitary wave prop- 
agates towards the plasma- vacuum interface against the 
plasma density gradient. Similarly to the 2D case dis- 
cussed in Refs. [hi], |l6"| , when the soliton approaches 
the plasma vacuum interface it radiates its trapped elec- 
tromagnetic wave away. This result shows clearly the 
difference between the wake field and the soliton. 

In conclusion we have demonstrated the existence of 
the three dimensional sub-cycle relativistic electromag- 
netic solitons in a collisionless cold plasma. The soli- 
tons consist of the electromagnetic and electrostatic fields 
with the structure of the oscillating electric dipole with 
the poloidal electric field and the toroidal magnetic field 
confined in the proplate cavity of the electron density. A 
substantial part of the pulse energy is transformed into 
solitons, approximately 25 — 30% of the incident laser 
pulse. The core of the soliton is positively charged on 
average in time and the soliton undergoes a Coulomb 
explosion in an ion time-scale. This process results in 
heating of the plasma ions. 
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